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Complementary Waveforms for Sidelobe Suppression and Radar Polarimetry 

Yuejie Chi, AH Pezeshki, and A. Robert Calderbank 


1. Introduction 

Advances in active sensing are enabled by the ability to control new degrees of 
freedom and each new generation of radar platforms requires fundamental advances 
in radar signal processing [1],[2], The advent of phased array radars and space-time 
adaptive processing have given radar designers the ability to make radars adaptable 
on receive [2], Modern radars are increasingly being equipped with arbitrary wave¬ 
form generators which enable transmission of different waveforms across space, time, 
frequency, and polarization on a pulse-by-pulse basis. The design of waveforms that 
effectively utilize the degrees of freedom available to current and future generation of 
radar systems is of fundamental importance. The complexity of the design problem 
motivates assembly of full waveforms from a number of components with smaller 
time-bandwidth products and complementary properties. By choosing to separate 
waveforms in space, time, frequency, polarization, or a combination of these degrees 
of freedom, we can modularize the design problem. Modularity allows us to first an¬ 
alyze the effect of each control variable in isolation and then integrate across degrees 
of freedom (e.g., time and polarization) to enable adaptive control of the radar’s oper¬ 
ation. 

Sixty years ago, efforts by Marcel Golay to improve the sensitivity of far infrared 
spectrometry led to the discovery of complementary sequences [3]—[5], which have 
the property that the sum of their autocorrelation functions vanishes at all delays other 
than zero. Almost a decade after their invention, Welti rediscovered complementary 
sequences (there are his D-codes) [6] and proposed to use them for pulsed radar. How¬ 
ever, since then they have found very limited application in radar as it soon became 
evident that the perfect autocorrelation property of complementary sequences cannot 
be easily utilized in practice. The reason, to quote Ducoff and Tietjen [7], is “in a 
practical application, the two sequences must be separated in time, frequency, or po¬ 
larization, which results in decorrelation of radar returns so that complete sidelobe 
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cancellation may not occur. Hence they have not been widely used in pulse compres¬ 
sion radars.” Various generalizations of complementary sequences including multiple 
complementary codes by Tseng and Liu [8] and multiphase (or polyphase) comple¬ 
mentary sequences by Sivaswami [9] and Frank [10] suffer from the same problem. 

In this chapter, we present a couple of examples to demonstrate how degrees of 
freedom available for radar transmission can be exploited, in conjunction with the 
complementary property of Golay sequences, to dramatically improve target detection 
performance. The chapter is organized as follows. In Section 2, we describe Golay 
complementary waveforms and highlight the major barrier in employing them for radar 
pulse compression. In Section 3, we will consider coordinating the transmission of 
Golay complementary waveforms in time (over pulse repetition intervals) to construct 
Doppler resilient pulse trains, for which the ambiguity function is effectively free of 
range sidelobes inside a desired Doppler interval. This construction demonstrates the 
value of waveform agility in time for achieving radar ambiguity responses not possible 
with conventional pulse trains. The materials presented in this section are based on 
work reported in [11]—[13]. In Section 4, we will present a new radar primitive that 
enables instantaneous radarpolarimetry at essentially no increase in signal processing 
complexity. In this primitive the transmission of Golay complementary waveforms 
is coordinated across orthogonal polarization channels and over time in an Alamouti 
fashion to construct a unitary polarization-time waveform matrix. The unitary property 
enables simple receive signal processing which makes the full polarization scattering 
matrix of a target available for detection on a pulse-by-pulse basis (instantaneously). 
The materials presented in this section are based on work reported in [14],[15],[11], 

2. Golay Complementary Waveforms 

2.1. Complementary Sequences 

Definition 1: Two length L unimodular 1 sequences of complex numbers x(t) and 
y{€) are Golay complementary if for k = — (L — 1), ..., (L — 1) the sum of their 
autocorrelation functions satisfies 

C x (k) + C y (k) = 2L5(k) ( 1 ) 

where C x (k) is the autocorrelation of x(£) at lag k and S(k) is the Kronecker delta 
function. This property is often expressed as 

\X(z)\ 2 + \Y(z)\ 2 = 2L (2) 

where X(z) is the z-transform of x(£) given by the polynomial 

L -1 

X(z) = J2 x(£)z~ e . (3) 

£=0 

1 The original complementary sequences invented by Golay for spectrometry were binary and polyphase 
complementary sequences were invented later. However, in this chapter we will not distinguish between 
them and refer to all complementary pairs as Golay complementary. 



Henceforth we may drop the discrete time index i from x(£) and y(£) and simply use 
x and y. Each member of the pair (. x , y) is called a Golay sequence. 

At first, it may appear that complementary sequences are ideal for radar pulse com¬ 
pression. However, as we will show in Section 2.2. the complementary property of 
these sequences is extremely sensitive to Doppler effect. For this reason, complemen¬ 
tary sequences have not found wide applicability in radar, despite being discovered 
almost sixty years ago. Some of the early articles that explore the use of Golay comple¬ 
mentary sequences (or their generalizations to polyphase and multiple complementary 
sequences) for radar are [6],[8]-[10], 

2.2. Golay Complementary Waveforms for Radar 

A baseband waveform w x (t) phase coded by a Golay sequence x is given by 

L—l 

w x (t ) = Y x(£)Q(t - £T C ) (4) 

where Q(t) is a unit energy pulse shape of duration T c and T c is the chip length. The 
ambiguity function \w x ("E v ) °f w x (t) is given by 

OO 

Xw x {t,v) = J w x (t)w x (t - T)e~ 3Vt dt 

— OO 

L—l 

= Y A x (k,uT c )xn(T + kT c ,u) (5) 

k=-(L- 1) 


where w(t) is the complex conjugate of w(t), xn(r, v) is the ambiguity function of 
the pulse shape f1(f), and A x (k, vT c ) is given by 

L — l 

A x (k, vT c ) = Y, x(l)x(£ - k)e~ 3 ^. (6) 

e=-(L- 1) 

If complementary waveforms w x (t) and w y (t) are transmitted separately in time, 
with a T sec time interval between the two transmissions, then the ambiguity function 
of the radar waveform w(t) = w x (t) + w y (t — T) is given by 2 

Xw(r, v) = Xw x (t, v) + e Jl ' T Xw y (r, v) ■ (7) 

Since the chip length T c is typically very small, the relative Doppler shift over chip 
intervals is negligible compared to the relative Doppler shift over T, and the ambiguity 

2 The ambiguity function of w(t) has two range aliases (cross terms) which are offset from the zero-delay 
axis by d zT. In this chapter, we ignore the range aliasing effects and refer to \w (t, v) as the ambiguity 
function of w(t). Range aliasing effects can be accounted for using standard techniques devised for this 
purpose (e.g. see [16]) and hence will not be further discussed. 
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Figure 1 . The perfect autocorrelation property of a Golay pair of phase coded wave¬ 
forms 

function \w (t, v) can be approximated by 

L—l 

Xw(t,v)= ^2 [C x {k) + e 3vT Cy{k)]xn{T+ kT c ,v) (8) 

k=-(L- 1 ) 

where in this approximation we have replaced A x (k 1 uT c ) and A y (k,vT c ) with the 
autocorrelation functions C x (k) and C y (k), respectively. 

Along the zero-Doppler axis (v = 0), the ambiguity function Xw( T ,v) reduces to 
the autocorrelation sum 

L—l 

Xw(t, 0) = ^2 [Cx{k) + Cy(k)\xn(r + kT c ,0) = 2L X a(r,0) (9) 

k— — (L—l) 

where the second equality follows from (1). This means that along the zero-Doppler 
axis the ambiguity function Xjh(t. v) is “free” of range sidelobes. 3 This perfect au¬ 
tocorrelation property (ambiguity response along zero-Doppler axis) is illustrated in 
Fig. 1 for a Golay pair of phase coded waveforms. 

Off the zero-Doppler axis however, the ambiguity function Xw(t, v) has large side- 
lobes in delay (range) as Fig. 2 shows. The color bar values are in dB. The range side- 

3 The shape of the autocorrelation function depends on the autocorrelation function \o (r. 0) for the 
pulse shape f2(f). The Golay complementary property eliminates range sidelobes caused by replicas of 
XQ (t, 0) at nonzero integer delays. 











Doppler(rad) 


Figure 2 . The ambiguity function of a Golay pair of phase coded waveforms separated 
in time 


lobes persist even when a pulse train is constructed in which w(t) = w x (t) + w y (t — T) 
is transmitted several times, with a pulse repetition interval (PRI) of T seconds be¬ 
tween consecutive transmissions. In other words, the ambiguity function of a conven¬ 
tional pulse train of Golay complementary waveforms, where the transmitter alternates 
between w x (t) and w y (t) during several PRIs, also has large range sidelobes along 
nonzero Dopplers. The presence of Doppler-induced range sidelobes in the ambiguity 
function means that a strong reflector can mask a nearby (in range) weak target that is 
moving at a different speed. Figure 3 shows the delay-Doppler map at the output of a 
radar receiver (matched filter), when a conventional (alternating) pulse train of Golay 
complementary waveforms was transmitted over 256 PRIs. The horizontal axis shows 
Doppler and the vertical axis depicts delay. Color bar values are in dB. The radar 
scene contains three stationary reflectors at different ranges and two slow-moving tar¬ 
gets, which are 30dB weaker than the stationary reflectors. We notice that the range 
sidelobes from the strong stationary reflectors make it difficult to resolve the weak 
slow-moving targets. The sensitivity of Golay complementary waveforms to Doppler 
effect has been the main barrier in adopting these waveforms for radar pulse compres¬ 
sion. As we will show in Section 3, this issue can be resolved by properly coordinating 
the transmission of Golay complementary waveforms over time. 

3. Doppler Resilient Pulse Trains 

Waveform agile radars are capable of changing their transmit waveform on a pulse- 
by-pulse basis, and hence are not limited to transmitting the same waveform (or wave- 





















































N=256, Alternating sequence with output in dB 
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Figure 3 . Delay-Doppler map at the output of a radar receiver (matched filter), when 
a conventional (alternating) pulse train ofGolay complementary waveforms is trans¬ 
mitted. 

form pair) over time. It is natural to ask whether or not it is possible to exploit this new 
degree of freedom to construct a Doppler resilient pulse train (or sequence) of Golay 
complementary waveforms, for which the range sidelobes of the pulse train ambigu¬ 
ity function vanish inside a desired Doppler interval. Following the developments in 
[11]—[13], we show in this section that this is indeed possible and can be achieved by 
coordinating the transmission of a Golay pair of phase coded waveforms ( w x (i),w y (t )) 
in time according to l’s and —l’s in a properly designed biphase sequence. 

Definition 2: Consider a biphase sequence V = { p n }n=o > Vn G { — 1,1} of length 
N, where N is even. Let 1 represent w x (t) and let —1 represent w y (t). The V-pulse 
train w-p(t) of (w x (t),w y (t)) is defined as 

1 JV_1 

w v (t) = - Y [(1 + p„)w x (t - nT) + (1 - p n )wy(t - nT )] . (10) 

n—0 

The nth entry in w-p(t) is w x (t) if p n = 1 and is w y (t) if p n = — 1. Consecutive 
entries in the pulse train are separated in time by a PRI T. 

The ambiguity function of the V -pulse train w-p(t), after ignoring the pulse shape 
ambiguity function and discretizing in delay (at chip intervals), can be written as [13] 

N-l N -1 

Xw v (k, 0) = ~[C x (k ) + Cy{k)] Y, e 3ne + iy[G x {k) - C y (k)} Y Pne jn6 ( 11 ) 

n—0 n—0 

























where 9 = vT is the relative Doppler shift over a PRI. The first term on the right-hand- 
side of (11) is free of range sidelobes due to the complementary property of Golay 
sequences x and y. The second term represents the range sidelobes as C x (k) — C y (k ) 
is not an impulse. We notice that the magnitude of the range sidelobes is proportional 
to the magnitude of the spectrum S-p(9) of the sequence V , which is given by 


N-l 

S V (d)=Y J Pn ein9 ■ ( 12 ) 

71=0 

By shaping the spectrum S-p(O) we can control the range sidelobes. The question is 
how to design the sequence V to suppress the range sidelobes along a desired Doppler 
interval. One way to accomplish this is to design V so that its spectrum S-p(9) has a 
high-order null at a Doppler frequency inside the desired interval. 

3.1. Resilience to Modest Doppler Shifts 

Let us first consider the design of a biphase sequence whose spectrum has a high- 
order null at zero Doppler frequency. The Taylor expansion of S-p (6) around 0 = 0 is 
given by 

OO 

= d3) 

i=0 n ' 

where the coefficients 0) are given by 


N-l 

/f(0) = i‘E4, * = 0,1,2,... (14) 

71=0 

To generate an Mth-order null at 9 = 0, we need to zero-force all derivatives 0) 

up to order M, that is, we need to design the biphase sequence V such that 

N-l 

J2n t Pn = 0, I-- 0.1.2...,,.!/. (15) 

71=0 

This is the famous Prouhet (or Prouhet-Tarry-Escott) problem [17],[18], which we 
now describe. 

Prouhet Problem: [17],[18] Let S = {0,1 ,N — 1} be the set of all integers 
between 0 and N— 1. Given an integer M, is it possible to partition S into two disjoint 
subsets So and Si such that 

E r m = E 7" 1 < 16 > 

t*£§o i 

for all 0 < to < All Prouhet proved that this is possible only when N = 2 M+1 and 
that the partitions are identified by the Prouhet-Thue-Morse sequence. 

Definition 3: [18],[19] The Prouhet-Thue-Morse (PTM) sequence V = (jpk)k >o 
over {—1,1} is defined by the following recursions: 



1- Po = 1 
2. P2k = Pk 
3- P2k+1 = —Pk 

for all k > 0. 

Theorem 1 (Prouhet) [17],[18]: Let V = {pk)k> o be the PTM sequence. Define 
S 0 = {r e S = {0,1,2,..., 2 M+1 — 1}| p r = 0} 

(17) 

Sr ={geS={0,l,2,... ) 2 M+1 - 1}| p q = 1} 

Then, (16) holds for all m, 0 < m < M. 

We have the following theorem. 

Theorem 2: Let V = {p n }r,-o^ be the length-2 M PTM sequence. Then, the spec¬ 
trum S-p(6) of V has an (M — l)th-order null at 9 = 0. Consequently, the ambiguity 
function of a PTM pulse train of Golay complementary waveforms (transmitted over 
2 m PRIs) has an (M — l)th-order null along the zero Doppler axis for all nonzero 
delays. 

Example 1: The PTM sequence of length 8 is 

V = (p fc )Lo = +l-l-l+l-l+l+l-l. (18) 

The corresponding pulse train of Golay complementary waveforms is given by 

wpTM{t) = w x (t) +w y (t-T) + w y (t - 2 T) + w x (t- 3 T) 

+w y (t — 4 T) + w x (t — 5 T) + w x (t — 6 T) + w y (t — 7T). (19) 

The ambiguity function of 'WpTM(t) has a second-order null along the zero-Doppler 
axis for all nonzero delays. 

Figure 4(a) shows the ambiguity function of a length -(N = 2 8 ) PTM pulse train 
of Golay complementary waveforms, which has a seventh-order null along the zero- 
Doppler axis. The horizonal axis is Doppler shift in rad and the vertical axis is delay in 
sec. The magnitude of the pulse train ambiguity function is color coded and presented 
in dB scale. A zoom-in around zero-Doppler is shown in Fig. 4(b). We notice that the 
range sidelobes inside the Doppler interval [—0.1,0.1]rad have been cleared out. They 
are at least 80 dB below the peak of the ambiguity function. This means that a PTM 
pulse train can resolve a weak target that is situated in range near a strong reflector, 
provided that the difference in speed between the weak and strong scatterers is rela¬ 
tively small. Figure 5 shows the ability of the PTM pulse train in bringing out weak 
slow-moving targets (e.g. pedestrians) in the presence of strong stationary reflectors 
(e.g. buildings) for the five target scenario discussed earlier in Fig. 3. This example 
clearly demonstrates the value of waveform agility over time for radar imaging. We 
note that the Golay complementary sequences used for phase coding in obtaining Figs. 



N=256, PTM with output in dB 



Figure 4 . Ambiguity function of a lengthfN = 2 8 ) PTM pulse train ofGolay comple¬ 
mentary waveforms: (a) the entire Doppler band (b) Doppler band [—0.1, 0.1 }rad 





































































































































N=256, PTM with output in dB 
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Figure 5 . The PTM pulse train of Golay complementary waveforms can bring out 
weak targets which would have otherwise been masked by the range sidelobes of 
nearby strong reflectors that move at slightly different speeds. 


3-5 are of size 64 and the pulse shape is a raised cosine. The chip length is T c = 100 
nsec, the carrier frequency is 17 GHz, and the PRI is T = 50 //sec. 

3.2. Resilience to Higher Doppler Frequencies 


We now consider the design of biphase sequences whose spectra have high-order 
nulls at Doppler shifts other than zero. Consider the Taylor expansion of S-p ( 8 ) around 
0 = 9 0 ± 0 : 

OO 1 

Sv(0) = y / ^f% ) (0 o )(e-6 o ) t ( 20 ) 

' n\ 

t =0 

where the coefficients /p \o 0 ) are given by 

JV-l 

= j t J2 nt P” ejn(>0 > * = 0 , 1 , 2 ,... ( 21 ) 

9—9q n —0 


$\0o) = 


df 


We wish to zero-force all the derivatives fp\6o) up to order M, that is, we wish to 
design the sequence V so that 

$(0o) = 0, for all t = 0,1,..., M. (22) 

Here, we limit our design to rational Doppler shifts do = 2 ttII ra , where l and m 1 
are co-prime integers and assume that the length of V is N = mq for some integer q. 













If we express 0 < n < TV— 1 as n = rm+i, where 0 < r < q— 1 and 0 < i < m— 1, 
then using binomial expansion for n 4 = (rm + i ) 4 we can write /p (0q) as 


/£ } (0o) 


N-l 


Jn0 o 


f ^2 rfpne 

n —0 

9 — 1 m— 1 

i‘EE {rm. + iYprm+ie 

r —0 i=0 
9—1 m—1 t 


t j{rm+i) 


S („)** “(rm)V m+a e J ' 

r—0 i=0 w—0 E / 

4 /a m_1 

Q m “ E j 

u—0 ' ' i= 0 


*-v" 


9-1 

^ Prm+i 
r —0 


(23) 

(24) 

(25) 

(26) 


Define a length-g sequence {fr r jY=o as b r = Prm+i, 0 < r < q — 1. If {b r }\Lj 
satisfies 

9-1 

^ r u b r = 0, for all 0 < u < t (27) 

r=0 


then the coefficient /p (0o) will be zero. From Theorem 1, it follows that the zero¬ 
forcing condition in (27) will be satisfied if {& r ]v=o * s the PTM sequence of length 2*. 
We note that o) is automatically zero as e= 0. Therefore, to zero- 

force the derivatives f^{0o) for all f < M, it is sufficient to select V = {p n }n= cT _1 
such that each i = 0, • • • , to— 1 is the length-2 M PTM sequence. We call 

such a sequence a (2 M ,m)-PTM sequence. The (2 m ,to)-PTM sequence has length 
2 m x to and is constructed from the length-2 M PTM sequence by repeating each 1 
and —1 in the PTM sequence m times, that is, it is constructed by oversampling the 
PTM sequence by a factor to. 


We have the following theorem. 

Theorem 3: Let V = {p n }n= cT _1 t> e th e (2 M , ?n)-PTM sequence, that is to say 
each of {prm+i}r=o~ 1 , * = 0, • • • , to — 1 is a PTM sequence of length 2 M . Then the 
spectrum S-p{9) of V has Mth-order nulls at all 9q = 2irl/m where l and to ^ 1 are 
co-prime integers. 

Corollary 1: Let V be the (2 m ,to)-PTM sequence. Then the spectrum S-p{9) of 
V has an (M — l)th-order null at 6*o = 0 and (M — h — l)th-order nulls at all 0q = 
2nl/{2 h m). where l and to ^ 1 are co-prime, and 1 < h < M — 1. 

Proof: The proof for 9q = 0 is straightforward. For 9 0 = 2irl/(2 h m), we have 


$\0o) = t 

u —0 


/.\ 2 h m— 1 

( j(2 h m) u ^2 

r4-_i "1 

^ ^ T P2 h mr+i 

4=0 

r= 0 


(28) 



where q = 2 M . The corollary follows from the fact that downsampling a PTM se¬ 
quence by a power of 2 produces a PTM sequence of shorter length. 

Theorem 3 and Corollary 1 imply that oversampled PTM pulse trains can bring out 
weak targets which are situated in range near strong reflectors, provided that the rela¬ 
tive distance in Doppler shift between the weak and strong scatterers is approximately 
2nl/(2 h m), where l and m ^ 1 are co-prime, and 1 < h <■ M — 1. 

Example 2: The spectrum of the (2 3 ,2)-PTM sequence 

7 ? = {+l,+l,—1,-1,—1,—1,+1,+1,—1,-1,+1,+1,+1,+1,—1,-1} 

has a third-order null at d 0 = n rad, a second-order null at 0 (t 0 rad, first-order nulls 

at 9q = 7t/ 2 rad and 9q = 3tt/2 rad, and zeroth-order nulls at do = (2 k + 1)tv/4 rad 
for k = 0,1,2,3. This spectrum is shown in Fig. 6 in solid black line. The spectrum 
of the (2 2 ,2)-PTM sequence is also shown in this figure in dashed black line. This 
spectrum has a second-order null at do = tt rad, a first-order null at 0q 0 rad, and 
zeroth-order nulls at do = 7r/2 rad and d 0 = 3tt/ 2 rad. 



Figure 6 . The spectra of ( 2 3 , 2)- and (2 , 2)-PTM sequences. 


Figure 7(a) shows the ambiguity function of a (2 s , 3)-PTM sequence of Golay com¬ 
plementary waveforms. The color bar values are in dB. This ambiguity function has 
an eighth-order null at d 0 = ±27t/ 3, a seventh-order null at zero Doppler, sixth-order 
nulls at do = ±7t/ 3, and so on. A zoom in around do = 27r/3 is provided in Fig. 7(b) 
to demonstrate that range sidelobes in this Doppler region are significantly suppressed. 
The range sidelobes in this region are at least 80 dB below the peak of the ambiguity 
function. This example shows the value of properly coordinating the transmission of 
complementary waveforms over time for achieving a desired ambiguity response. 
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Figure 7. Ambiguity function of the (2 8 , 3 )-PTM pulse train of Golay complementary 
waveforms: (a) the entire Doppler band (b) Zoom in around 9 q = 2n/3. 








































































































































































4. Instantaneous Radar Polarimetry 


Fully polarimetric radar systems are capable of transmitting and receiving on two 
orthogonal polarizations simultaneously. The combined signal then has an electric 
field vector that is modulated both in direction and amplitude by the waveforms on the 
two polarization channels, and the receiver is used to obtain both polarization compo¬ 
nents of the reflected waveform. The use of two orthogonal polarizations increases the 
degrees of freedom available for target detection and can result in significant improve¬ 
ment in detection performance. 

The radar cross section of an extended target such as an aircraft or a ship is highly 
sensitive to the angle of incidence and angle of view of the sensor (see [20] Sections 
2.7-8). In general the reflection properties that apply to each polarization component 
are also different, and indeed reflection can change the direction of polarization. Thus, 
polarimetric radars are able to obtain the scattering tensor of a target 


S = ( avv , (29) 

\<jhv &hhJ 

where <jvh denotes the target scattering coefficient into the vertical polarization chan¬ 
nel due to a horizontally polarized incident field. Target detection is enhanced by con¬ 
current rather than serial access to the cross-polarization components of the scattering 
tensor, which varies more rapidly in standard radar models used in target detection 
and tracking [21] than in models used in remote sensing or synthetic aperture radar 
[22],[23]. In fact what is measured is the combination of three matrices 

n=( hvv hvH \ = C Rx SC Tx , (30) 

\n HV riHH J 


where Cr x and Ct x correspond to the polarization coupling properties of the trans¬ 
mit and receive antennas, whereas S results from the target. In most radar systems 
the transmit and receive antennas are common, and so the matrices Ct x and Cr x are 
conjugate. The cross-coupling terms in the antenna polarization matrices are clearly 
frequency and antenna geometry dependent but for the linearly polarized case this 
value is typically no better than about -20dB. 

In this section, we follow [14],[15] to describe a new approach to radar polarime¬ 
try that uses orthogonal polarization modes to provide essentially independent chan¬ 
nels for viewing a target, and achieve diversity gain. Unlike conventional radar po¬ 
larimetry, where polarized waveforms are transmitted sequentially and processed non- 
coherently, the approach in [14],[15] allows for instantaneous radar polarimetry (IRP), 
where polarization modes are combined coherently on a pulse-by-pulse basis. Instan¬ 
taneous radar polarimetry enables detection based on full polarimetric properties of 
the target and hence can provide better discrimination against clutter. When compared 
to a radar system with a singly-polarized transmitter and a singly-polarized receiver 
instantaneous radar polarimetry can achieve the same detection performance (same 
false alarm and detection probabilities) with a substantially smaller transmit energy, or 
alternatively it can detect at substantially greater ranges for a given transmit energy. In 



this section, we present the main idea for waveform transmission in IRP. The reader is 
referred to [14],[15] for simulation studies that demonstrate the value of IRP. 


4.1. Unitary Waveform Matrices 

Let us employ both polarization modes to transmit four phase-coded waveforms wjj, 
Wy, w 2 H , Wy. On each polarization mode we transmit two waveforms separated by a 
PRI T. We employ Alamouti coding [24] to coordinate the transmission of waveforms 
over the vertical (V) and horizontal (H) polarization channels; that is, we define 

w 2 h = wh (31) 

Wy = -wjj, (32) 

where r denotes complex conjugate time-reversal. After discretizing (at chip intervals) 
and converting time-indexed sequences to ^-transform domain, we can write the radar 
receive equation in matrix form as 

R(z) = + Z (z) (33) 


where RO) is the 2 by 2 radar measurement matrix at the receiver, H is the 2 by 2 
scattering matrix in (30) and Z (z) is a noise matrix. The entries of H are taken to be 
constant since they correspond to a fixed range (delay d ) and a fixed time. For now, 
we limit our analysis to zero-Doppler axis and will postpone the treatment of Doppler 
effect to Section 4.2. The Alamouti waveform matrix W (z) is given by 


wo) 


(wfo) -Whiz)] 
\Whiz) W^z) ) 


(34) 


where W(z) = z L Wiz x ) for a length L sequence w. 


If we require the matrix W 0) to be unitary, that is. 


W(z)W(z) 


Whiz) -Whiz)] ( Whiz) Whiz)\ =2L L (1 0 
[Wh(z) Whiz) ) [-Whiz) Whiz)) ~ [o 1 

(35) 


then it is easy to estimate the scattering matrix H by post-multiplying (33) with W (z). 
The unitary condition is equivalent to 


\Whiz)\ 2 +\Whiz)\ 2 = 2L. 


(36) 


This is the same condition as the Golay complementary condition and is satisfied 
by choosing Wy(z) = X(z) and WOO) = Y0)> where Xiz) and YO) are 21- 
transforms of the Golay complementary sequences x and y, respectively. This shows 
that by properly coordinating the transmission of Golay complementary waveforms 
across polarizations and over time we can make the four channels HH, VV, VH, and 
HV available at the receiver (with delay L) using only linear processing. The four 
matched filters (in z-domain) are given by 

ntA - ( Ml ( z )Rv(z) M 2 (z)R v (z)\ 

^ Z) ~ [M^Rniz) M 2 iz)R H iz)J 


(37) 



where 


M 1 (z) = W^z)z~ L -W^{z) 

M 2 (z) = Wh(z)z~ L + Wy(z). (38) 

and Rv{z) denotes the first row of li.(z) and represents the 2 -transform of the entire 
radar return (over two time slots) received on the vertical polarization channel. 


Remark 1: The above description suggests that a radar image will be available only 
on every second pulse, since two PRIs are required to form an image. However, after 
the transmission of the first pulse, images can be made available at every PRI. This is 
done by reversing the roles of the waveforms transmitted on the two pulses. Thus, in 
the analysis, the matrix W(z) in (33) is replaced by 


V(2) 


(-W^(z) W^)\ 

\W^(z) Wh(z)J 


(39) 


which is still unitary due to the interplay between Golay property and Alamouti cod¬ 
ing. Moreover, the processing involved is essentially invariant from pulse to pulse: the 
return pulse in each of the and channels is correlated against the transmit pulse on that 
channel. This yields an estimate of the scattering matrix on each pulse. 


4.2. Doppler Resilient IRP 

In the previous section, we restricted our analysis to zero-Doppler axis. Off the zero- 
Doppler axis, a relative Doppler shift of 9 exists between consecutive waveforms on 
each polarization channel and the radar measurement equation changes to 

RO) = r d HW(z)D(0) + Z(z) (40) 

where D(6>) = diag(l, e- ?e ). Consequently, the unitary property of W(z) can no 
longer be used to estimate H, due to the fact that W( 2 )D(#)W( 2 ) is not a factor of 
identity. In other words, off the zero-Doppler axis the four elements in W (i:)D(0) W(z) 
have range sidelobes. 

The matrix W(2)D(0)W(2) can be viewed as a matrix-valued ambiguity function 
(in z-domain ) for IRP. In Section 3, we showed that the range sidelobes of the am¬ 
biguity function of a pulse train of Golay waveforms can be suppressed (in a desired 
Doppler band) by carefully selecting the order in which the Golay complementary 
waveforms are transmitted. This result can be extended to IRP. In fact, it is shown 
in [11] that the range sidelobes of the IRP matrix-valued ambiguity function can be 
cleared out along modest Doppler shifts by coordinating the transmission of the fol¬ 
lowing Alamouti matrices 


Xi(z) 


(X(z) -J(*)\ 

1^(2) X(z) ) 


and X_i ( 2 ) 


(-Y(z) -X(z)\ 

V *(*) -nz)) 


(41) 


according to the Is and —Is in the PTM sequence, where 1 represents Xi( 2 ) and 
— 1 represents X_i( 2 ), and X(z) and Y ( 2 ) are 2 -transforms of the Golay pair (a;, y). 



Range sidelobes along higher Doppler intervals can be cleared out by coordinating the 
transmission of Xi(^) and X (z) according to oversampled PTM pulse trains. 

5. Conclusion 

In this chapter, we focused on the use and control of degrees of freedom in the radar 
illumination pattern and presented examples to highlight the value of properly utilizing 
degrees of freedom. We showed that by coordinating the transmission of Golay com¬ 
plementary waveforms in time (or exploiting waveform agility over time) according to 
carefully designed biphase sequences we can construct pulse trains whose ambiguity 
functions have desired properties. We also showed that by combining Alamouti coding 
and Golay complementary property we can construct unitary polarization-time wave¬ 
form matrices that make the full polarization scattering matrix of a target available for 
detection on a pulse-by-pulse basis. Looking to the future, we see unitary waveform 
matrices as a new illumination paradigm that enables broad waveform adaptability 
across time, space, frequency and polarization. 
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